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I.  INTRODUCTION 


In  hydrography,  the  character  of  the  bottoin  should  be  determined  for  nautical 
charting,  particularly  in  harbors,  anchorage  areas  or  other  areas  where  the  bottom 
characteristic  is  a  significant  factor  to  be  considered  for  navigation  safety. 

Conventional  sampling  methods,  using  clamshell  snappers  (L’mbach,  1976).  are  used  to 
define  the  characteristics  of  the  bottom  surface  layer.  It  is  tedious  and  time 
consuming,  especially  in  off-shore  geologic  investigations.  If  some  kind  of  theory  or 
model  existed  to  remotely  classify  the  sediments,  the  task  could  be  carried  out  easier 
and  faster.  Moreover,  this  ability  would  permit  nautical  charts  to  contain  much  more 
information. 

In  the  ocean  bottom,  the  magnitude  and  nature  of  the  sound  scattering  are 
functions  of  both  the  particle  size  of  the  bottom  and  surface  (bottom)  relief  (Mckinney 
and  Anderson.  296-).  The  sound  echo  from  sediments  is  subject  to  three  types  of 
reflecting  or  scattering  ^Clark.  Proni.  Seem,  Tsai.  19S5)  : 

1)  coherent  reflection  from  the  water-sediment  interface;  the  signal  amplitude 
will  be  approximately  that  given  by  spherical  spreading  of  the  wave  front 
from  the  image  source  modified  by  the  plane  wave  reflection  coefficient 
(Clay,  1966), 

2)  incoherent  or  statistically  variable  scattering  from  the  surface  irregularities 
of  the  sediment,  and 


3)  incoherent  scattering  from  within  the  volume  of  the  sediment  caused  by  the 
acoustically  irregular  matrix  of  the  sediment. 

Items  1  and  2  arise  from  the  surface  of  the  sediments  and  item  3  from  the  sound 
pulse  penetrating  into  the  sediment  volume  and  then  being  scattered  back  to  the 
transducer  by  the  granular  matrix  of  the  sediment.  Since  different  sediments  have 
different  acoustic  matrices,  they  have  different  volume  reverberation  characteristics  and 
there  should  exist  a  mathematical  relationship  between  sediment  properties  and  the 
volume  reverberation. 

The  major  goal  of  this  experiment  was  to  establish  a  fundamental  laboratory 
procedure  for  the  remote  sensing  of  ocean  sediment  volume  reverberation.  This  was  to 
be  accomplished  by  accurately  measuring  the  properties  of  a  pulse  reflected  from  the 
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sediment.  In  application,  since  the  reflected  sound  puise  is  received  by  the  same 
transducer,  i.e.,  "monostatic  sonar",  our  discussion  is  restricted  to  "back-reflection"  and 
"backscattering". 
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II.  BACKGROUND 


Bradshaw  (1981)  used  Monterey  #31  fine  sand,  fne  same  material  used  in  this 
experiment,  to  study  the  sound  propagation  into  a  fast  bottom  medium.  In  his  study, 
air  bubbles  were  removed  from  within  the  sediment  based  on  heating  the  water  sand 
system.  Bleach  was  added  at  a  ratio  of  one-half  gallon  bleach  to  70  gallons  of  water 
plus  sand  in  order  to  control  biologic  growth.  The  physical  properties,  such  as  density, 
porosity,  speed  of  sound,  reflecton  coefficient  etc.,  were  measured.  Kosnik  (1984)  used 
a  high  pressure  water  jet  to  remove  the  air  bubbles  from  the  water,  sand  mixture. 
However,  no  volume  scattering  measurements  were  made  by  either  investigator. 

Morse  and  Ingard  (1968)  give  a  theoretical  derivation  on  incoherent  scattering. 
Under  the  assumptions  that  :  I)  the  scattering  objects  are  quite  small  compared  with 
wave  length  X.  2)  the  dimensions  of  scattering  region  are  much  larger  than  X,  3  )  the 
scatterers  are  all  spheres,  and  4)  the  scatterers  are  populated  sparsely  and  no  mutiple 
scattering  existed,  the  incoherently  scattered  intensity  1-  per  unit  incident  intensity  I  - , 
at  angle  8  to  the  incident  wave,  for  frequency  o)/2tt,  is 


(2.1) 


(J_ K"  ar  siiM(l  2)0)  exp(-2K“  a"  siiC(  1/2)0) 
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where  X  is  the  total  number  of  scatterers  within  unit  volume,  K  is  the  wave  number  for 
the  mixture,  a  is  the  radius  of  scatterers.  Furthermore,  for  monostatic  sonar,  the 
volume  scattering  coefficient  for  unit  volume  can  be  reduced  as 


(2.2) 


where 


K 


n 


K 


0 


K 


0 
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3  Pn  -  3  p0 

Yp  = 

2  Pn  *  Po 

and  Kn  is  the  compressibility  of  the  scattered  kq  is  the  compressibility  of  the  fluid.  pu 
is  the  density  of  the  scatterer  and  p_  is  the  density  of  fluid. 

For  a  small  (ka  <  <  1)  rigid  sphere,  the  backscattering  cross  section  is  given 
(Urick,  1983)  as 

(T  =  2.S  ( 7t  a2){K  a)4 

Note  that,  as  a  single  sphere,  it  must  be  a  single  scattering.  The  scattering  coefficient 
for  unit  volume  can  be  obtained  as 

sv  =  N  a  (2.3) 

Both  theories  are  not  very  appropriate  for  the  real  ocean  sediment  volume 
reverberation.  However,  these  are  the  only  two  theories  available.  The  backscattering 
of  sound  from  the  ocean  bottom  has' been  studied  for  many  years,  but  unfortunately  no 
papers  were  found  to  discuss  the  volume  reverberation  of  ocean  sediments. 
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III.  DESCRIPTION  OF  VOLUME  REVERBERATION 


As  sound  propagates  in  a  fluid,  part  of  the  sound  is  intercepted  and  reradiated  by 
inhomogeneites  in  the  fluid.  The  reradiation  of  sound  is  called  scattering.  The  total  of 
the  scattering  contributions  from  all  the  scattercrs  is  called  reverberation  (Crick.  1983). 

Volume  reverberation  is  defined  as  sound  scattered  back  to  the  transducer  by 
scattering  centers  in  the  volume  of  the  sea  (National  Defense  Research  Commitee, 
1969).  Volume  reverberation  may  be  considered  as  a  blending  of  large  number  of 
echoes.  It  has  been  studied  mostly  in  the  Deep  Scattering  layer  .  a  region  of  high 
biological  activity.  Reverberation  level,  in  decibels,  is  used  to  measure  volume 
reverberation.  Consider  a  directional  projector  in  an  ocean  containing  a  large  number 
of  volume  scatterers  (Figure  3.1).  Let  the  axial  intensity  at  unit  distance  be  IQ.  At  a 
distance  r  where  there  is  a  volume  element  dV  of  volume  scattercrs,  the  incident 
intensity  is 

=  J0  f2 

The  scattered  intensity  i  m  from  the  volume  dV  is 
dls(l>=  Ij  sv  dV 

where  s,.  is  the  volume  scattering  coefficient.  The  scattered  intensity  back  at  the 
receiver  is 

dl$  (r)  =  Ij  sv  dV.r2 

=  (I0.r4)svdV 

and  the  scattered  intensity  from  all  scatterers  contributing  to  the  scattered  intensity  at 
time  t  is 

Is(0  =  0o  A  sv  V 

where  V  is  the  total  volume  contributing  for  the  geometry  of  Figure  3.1  calculated  by 
(Crick,  19S3) 

V  =  .Qr2  ct/2 
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TRANSDUCER 


Figure  3.1  Portion  of  ocean  scattering. 


IS 


where  Cl  is  the  solid  angle  of  the  sound  beam,  x  is  the  pulse  length  and  c  is  the  sound 
speed,  so 


is  =  n0  f2)  sv  ^  CT  2 
The  reverberation  level  is 

RL  =  10  LOG  Is(r).Tref 

=  SL  -  TL  -  Sv  +  10  LOG  (O  ct  2) 


where  TL  is  the  transmission  loss  for  spherical  spreading  20LOG(r)  and  Sy  is  the 
scattering  strength  for  unit  volume  10LOG(sy).  In  Deep  Scattering  Layer  studies,  this 
reverberation  level  equation  is  sufficient,  because  the  received  pulse  is  regarded  as 
totally  coming  from  volume  reverberation,  i.e.  negligible  reflection  and  surface 
scattering.  Furthermore,  the  derivation  of  V  is  based  on  a  square  acoustic  pulse.  But 
for  a  real  transducer,  the  square  pulse  is  transformed  into  a  pulse  with  an  exponential 
rise,  a  flattened  top  and  an  exponential  decay.  The  sound  pulse  reflected  from  the 
sediment  and  received  by  the  transducer  combines  coherent  reflected,  incoherent 
scattering  from  the  surface  and  volume  reverberation. 

Let  a  transducer  located  in  the  water  send  out  a  sound  pulse  to  a  smooth 
water  air  interface.  The  received  echo  waveform  is  purely  reflection.  If  the  water 
interface  is  rough,  then  the  received  echo  waveform  includes  both  reflection  and  surface 
scattering.  Since  sound  energy  is  scattered  by  the  rough  surface,  the  received  waveform 
amplitude  is  reduced.  Consider  first-order  scattering  and  suppose  that  the  surface 
roughness  is  the  same  for  every  position.  The  amplitude  of  the  combined  sound  wave 
(reflected  and  scattered)  is  reduced  proportional  to  the  purely  reflected  wave  amplitude. 
In  Figure  3.2(a),  the  solid  line  represents  an  ensemble  of  reflected  waveform,  and  the 
dashed  line  represents  an  ensemble  of  reflected  and  scattered  waveform. 

The  contributions  to  a  pulse  reflected  from  a  sediment  surface  is  shown  on 
Figure  3.2(b).  The  dashed  line  is  the  contribution  by  reflection  and  surface  scattering: 
again,  only  first  order  scattering  is  considered.  The  dotted  line  includes  volume 
reverberation.  The  level  of  the  signal  is  raised  above  that  of  a  dashed  line  because  of 
energy  returned  from  within  the  sediment.  In  particular,  the  tail  decays  more  slowly 
than  for  reflection  from  the  water  air  interface.  This  is  the  significant  feature  of 
sediment  volume  reverberation  that  we  are  after  in  this  paper. 
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Figure  3.2  Theoretical  echo  waveform  for  water  and  sediment. 
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Two  alternative  assumptions  have  to  be  made  for  reflection  from  the  sediment 
surface  :  either  !)  the  surface  scattering  is  negligible  compared  with  volume 
reverberation,  or  2)  if  it  is  not  negligible,  it  must  have  the  same  amplitude  factor  as 
discussed  in  Figure  3.2(b). 

According  to  these  assumptions,  if  a  sediment  echo  is  subtracted  from  the  water 
echo  it  will  give  a  measure  of  the  volume  reverberation  in  the  sediment  (Figure  3.3). 

The  decaying  tail  of  the  echo  reflected  from  the  smooth  water  air  interface  is  an 
exponential  decay. 

vu.  (.)  -  v0  e'5t 

where  \’w  (t)  is  the  received  voltage  at  time  t,  5  is  the  decay  constant.  V  is  a  constant 
voltage,  and  t  is  the  time  starting  from  the  beginning  of  decay.  If  we  assume  the 
sediment  echo  also  decays  exponentially,  then  we  will  obtain  a  similar  function 

Vs  (t)  -  V0e-’l' 

where  V  (t)  is  the  received  voltage  from  sediment  at  time  t  and  i]  is  a  decay  constant. 
If  \’w  and  V$  are  both  normalized  to  the  same  initial  value  K.  the  difference  becomes 

Vv  (t)  =  K  (e'W  -  e‘5t)  (3.1) 


If  both  t’(t  and  Ot  are  much  less  than  unity, 

Vv(t)  =  K  (  n  -  6  )t  t  <  <  1  n  and  1  6 
This  initial  decay  is  linear  and  depends  on  the  difference  6-t|. 
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ECHO  FROM  WATER  SURFACE 


ECHO  FROM  SEDIMENT  SURFACE 


Figure  3.3  Real  echo  waveform  displayed  on  CRT. 
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IV.  EXPERIMENT  DESIGN 


A.  MATERIAL  AND  TANK  SELECTED 

Fresh  (tap)  water,  t?30  Monterey  fine  sand  and  Monterey  Aquarium  =2  for  a 
rough  gravel  sediment  were  the  media  used  in  the  experiment.  The  grain  size  of  =30 
Fine  sand  varies  from  0.15  mm  to  0.7  mm  (Bradshaw,  1981)  with  an  average  size  0.3 
mm.  The  Aquarium  =2  has  a  mean  size  5.3  mm  (Diaz.  19S6).  The  sand  water  mixture 
was  repeatedly  stirred  to  remove  air  bubble  trapped  while  transferring  the  sand  to 
different  tanks. 

Initial  experiments  were  conducted  in  the  NFS  aneehoic  tank  to  measure  the 
reflection  from  the  water/air  surface,  with  the  transducer  located  at  different  depth. 

The  second  tank,  a  wooden  tank  measuring  80  cm  x  80  cm  x  60  cm,  was  constructed  to 
hold  about  30  cm  of  aggregate.  This  tank  was  lowered  and  suspended  85  cm  below  the 
water  air  surface  of  the  aneehoic  tank.  A  third  tank,  a  steel  glass  tank,  measuring  70 
cm  x  70  cm  x  60  cm  was  filled  with  24  cm  of  fine  sand  or  aggregate  with  36  cm  of 
water  above. 

B.  ELECTRONIC  EQUIPMENT 

A  schematic  drawing  of  the  equipment  is  shown  in  Figure  4.1.  Except  a  high 
pass  filter,  all  components  were  off-the-shelf.  General  Radio  model  1310  oscillator 
with  a  frequency  of  182.0  kHz  was  fed  simultaneously  into  a  Hewlett-Packard  5233L 
frequency  counter,  and  a  General  Radio  Type  396-A  tone  burst  generator  was  used  to 
generate  a  16  cycle  pulse.  Signal  output  from  Hewlett-Packard  IIP  467-A  amplifier 
was  fed  through  a  Datasonic  Transmit  Receive  (T/R)  switch  and  a  small  21  to  2S  kHz 
high  pass  filter,  then  into  the  F-41  transducer.  The  F-41  transducer  was  discussed  in 
detail  by  both  Borchardt  (1985)  and  Diaz  (1986).  Signals  returned  from  the  sediment 
were  amplified  20  dB  by  a  Hewlett-Packard  465-A  pre-amplifier,  then  passed  through  a 
Spencer-Kenncdy  Laboratories,  Inc.  model  302  variable  electronic  filter  (set  at  135  kHz 
high  pass)  to  eliminate  low  frequency  noise  before  being  passed  to  the  digital  and 
analog  oscilloscopes. 

The  waveform  displayed  on  309 1  Nicolct  Digital  Oscilloscope  was  sampled  and 
sent  as  a  stream  of  characters  to  the  computer.  To  reduce  manual  operations  on  the 
Nicolct  3091  during  the  experiment,  a  1  IP-3421  A  Data  Acquisition  and  Control  Unit 
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Figure  4. 1  Schematic  of  equipment. 
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was  adopted  to  receive  commands  from  a  HP-S6  computer  and  then  automatically 
operate  the  Nicolet  3091.  However,  it  was  noticed  that  if  the  Nicolet  3091  was  not 
"warmed  up"  enough,  the  DC-offset  jumped  up  and  down,  and  continual  manual 
adjustment  was  needed.  After  2  or  3  hours  of  running,  it  became  stable  and  no  further 
adjustment  was  required. 


C.  VERIFICATION  OF  MEASUREMENTS 

No  matter  how  deep  the  transducer  is  located  below  the  water  air  surface,  as 
long  as  the  transducer  is  level,  the  exponential  decay  of  the  tail  should  be  the  same.  In 
Diaz's  thesis,  the  measurements  were  inconclusive  because  the  transducer  was  fixed  at 


one  position  throughout  the  experiment.  This  was  not  adequate  to  prove  the  system 
functioned  properly.  To  prove  that  the  previous  experiments  were  run  properly  and 
results  were  obtained  correctly,  a  verification  was  attempted  by  measuring  the 
reflection  from  the  water  air  interface  with  the  same  procedure,  but  with  the  transducer 
at  various  depths.  Data  set  20  (Figure  4.2)  for  tiie  transducer  70  cm  below  the  surface, 
gives  a  slope  of -3. ~5  x  10'“  Np  ’is.  Data  set  2!  (Figure  4.3)  was  obtained  by  raising 
the  transducer  15  cm  and  the  slope  is  -3.77  x  10'“  Np  Jis.  Both  results  were  close  to 
the  previous  results  obtained  by  Diaz  (1986).  After  these  two  data  sets  were  collected, 
a  new  transducer  mounting  was  built  and  adopted  to  improve  the  experimental 
precision  and  operation  ability.  The  transducer  was  leveled  by  a  counter-weight  on  the 
new  carriage.  To  ensure  the  system  and  operation  worked  properly,  more  locations 
were  tried  to  obtain  the  decay  constant  of  the  water/air  surface. 

Data  sets  23  and  24  were  obtained  at  the  same  transducer  depth,  but  seperated 
40  cm  horizontally.  Data  sets  24  and  25  were  obtained  at  the  same  horizontal  position 
but  with  50  cm  difference  in  depth.  Data  set  26  was  acquired  by  changing  10  cm  in 
both  the  depth  and  position  from  data  set  25.  Data  sets  26  to  2S  resulted  from 
changing  the  depth  from  70  cm  to  30  cm  below  the  water  air  interface.  The  results  of 
these  experiments  are  listed  in  Table  1  and  show  consistencies  between  each  other. 
Waveform  ensembles  of  the  tail  for  each  data  set  are  plotted  on  Figure  4.4  to  Figure 
4.9. 


The  average  slope  -3.77  ±  0.077  x  10'“  Np  jis  is  smaller  than  -3.S4  ±  0. 14  x 
10'“  Np  Jis  which  was  measured  in  Diaz's  thesis.  Statistically,  these  two  slopes  are  the 
same.  The  smaller  standard  deviation  might  suggest  that  the  new  carriage  has  achieved 
a  better  precision. 
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Figure  4.2  Decay  envelope  for  reflection  from  water, 'air  interface, 
transducer  70  cm  below  the  surface  (data  set  20). 
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Figure  4.3  Decay  envelope  for  rcllcction  from  water, 'air  interface, 
transducer  55  cm  below  the  surface  (data  set  21). 
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TABLE  1 

WATER  AIR  REFLEC  TION 


DATA 

SLOPE 

CORRELATION 

SET 

(xiO“2) 

COEFFICIENT 

DATA20 

-3.  75 

-0. 9984 

DATA2 1 

-3.  77 

-0. 9795 

DATA22 

-3.  86 

-0. 9979 

DATA23 

-3.  62 

-0.  9724 

DATA24 

-3.  68 

-0. 9895 

DATA25 

-3.  86 

-0. 9956 

DATA2 6 

-3.  80 

-0. 9679 

DATA2 7 

-3.  84 

-0.  9888 

DATA28 

-3.  79 

-0. 9810 

MEAN  SLOPE  = 

-3.  77  x  10“2  Np/fts 

<7  =  0.  07  x  10~2  Np/Jis 


After  dealing  with  the  water  air  surface,  the  transducer  was  turned  over  and 
pointed  downward  toward  the  sediment  which  was  contained  within  the  80  cm  x  SO  cm 
x  60  cm  wooden  tray.  Different  transducer  positions  were  choosen.  Data  sets  31  to  33 
were  at  the  same  depth,  but  horizontally  3  cm  away  from  each  other.  The  results  show 
the  waveform  were  totally  inconsistant.  Data  sets  33  to  35  were  at  the  same  horizontal 
position  but  different  depths;  results  were  also  inconsistant.  Above  results  are 
presented  on  Table  2  and  plotted  on  Figure  4.10  through  Figure  4.14.  So  far, 
experiments  reveal  the  lack,  of  isotropic  homogeniety  of  the  volume  contributes 
reverberation. 

D.  MEASUREMENT  CONSIDERATION 

During  the  course  of  these  experiments  some  other  effects  were  considered  and 
studied  to  ensure  that  the  experiment  was  working  properly. 
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Figure  4.4  Dccav  envelope  for  rellection  from  water/air  interface, 
transducer  was  30  cm  below  surface  (data  set  23). 
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Figure  4.5  Decay  envelope  for  reflection  from  water/air  interface, 
trnsducer  was  30  cm  below  surface  (data  set  24). 
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CORRELATION  C0EFF-  -.9895 


Figure  4.6  Decay  envelope  for  reflection  from  water, 'air  interface, 
transducer  was  Si)  cm  below  surface  (data  set  25). 
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Figure  4.7  Decay  envelope  for  reflection  from  water/air  interface, 
transducer  was  70  cm  below  surface  (data  set  26). 
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Figure  4.8  Decav  envelope  lor  reflection  from  watcr/air  interlace, 
transducer  was  50  cm  below  surface  (data  set  27). 
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Figure  4.9  Dccav  envelope  for  reflection  from  watcr/air  interface, 
transducer  was  3U  cm  below  surface  (data  set  28). 
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TABLE  2 

SEDIMENT  WATER  REFLECTION  AND  SCATTERING 


DATA 

SLOPE 

CORRELATION 

SET 

(x!0“2) 

COEFFICIENT 

DATA3 1 

-2.20 

-0.  9138 

DAT A3 2 

-2.  98 

-0.  8908 

DATA33 

03 

i — 1 

1 

-0.  7838 

DAT A3 4 

1 

CJ 

o 

03 

-0.  8174 

DATA35 

-2.  90 

-0.  9844 

MEAN  SLOPE  =  -2.62  ; 

t  10“2  Mp/)is 

cr  =  0.  85  x  10"2  Np/us 

1.  Effects  of  Sampling  Rate 

The  sound  signal  being  used  was  fixed  at  1S2  kHz  (5.5  jls  period).  The 
Nicolet  3091  digital  oscilloscope  has  as  its  fastest  sampling  rate  1  ps  which  is 
inadequate  to  give  a  precise  reproduction  of  the  numeric  waveform.  However, 
observations  have  shown  "jitter"  in  the  triggering  of  the  scope,  so  for  consecutive 
pulses  the  scope  would  sample  different  parts  of  the  waveform.  To  obtain  an  accurate 
waveform,  a  large  number  of  waveform  were  sampled  and  averaged.  If  sufficient 
"jitter"  is  present,  the  voltage  measured  in  a  given  bin  should  van;  between  the 
maximum  and  minimum  voltage  in  the  waveform  in  the  vicinity  of  the  corresponding 
time.  Since  the  sine  wave  varies  between  minus  and  plus,  the  voltage  values  were  • 
squared  before  averaging  in  order  to  give  an  accurate  representation  of  the  waveform 
envelope. 

To  determine  the  maximum  number  of  samples  required  for  each  data  set 
(average  waveform),  the  water  air  interface  was  measured  with  different  number  of 
samples.  The  results  showed  that  50  samples  for  each  data  set  were  enough  to  give  a 
true  waveform. 

2.  Effects  of  Bottom  Reflection  Interference 

For  a  16  cycle  signal  with  a  frequency  182  kHz,  the  pulse  length  is  88  )ls. 

Since  the  sound  pulse  emitted  had  a  decay  time  of  approximately  100  ns,  the  echo  from 
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Figure  4.10  Dccav  envelope  for  reflection  from  aggregate,1 ’water 
interlace,  transducer  was  30  cm  above  surface  (data  set  31). 
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Figure  4.1 1  Decay  envelope  for  reflection  from  aegregatc/water 
interface,  transducer  was  30  cm  above  surface  (data  set  32). 
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Figure  4.12  Decay  envelope  for  reflection  from  aggregate, 'water 
intei  (ace,  transducer  was  30  cm  above  surface  (dataset  33). 
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Figure  4.13  Dccav  envelope  for  reflection  from  aggregate, 'water 
"interlace,  transducer  was  50  cm  above  surface  (data  set  34). 
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Figure  4.14  Decay  envelope  for  reflection  from  aggregate, 'water 
interlace,  transducer  was  70  cm  above  surface  (data'set  35). 
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the  smooth  water  air  interface  should  also  have  the  same  decay  time.  Since  the 
purpose  of  this  project  was  to  compare  the  decay  portion  of  echos  reflected  from  water 
and  sediments,  the  waveform  was  sampled  for  a  100  fis  interval  from  the  beaming  of 
decay.  The  sound  speed  in  the  aggregate  was  measured  to  be  1555  m  s.  For  the  fine 
sand.  Bradshaw  (19S1)  obtained  the  speed  of  sound  to  be  1610  m  s.  Combining  these 
sound  speeds  and  the  time  interval  1SS  us,  the  depth  of  sediment  that  contrubutes  to 
the  main  echo  is  1-1.6  cm  for  aggregate  and  15.1  cm  for  fine  sand.  Interference  caused 
by  reflection  from  the  bottom  of  the  tank  will  not  affect  the  experiment  as  long  as  the 
sand  is  thicker  than  15.1  cm  and  the  aggregate  is  thicker  than  Id. 6  cm.  Throughout 
our  experiments,  we  kept  both  sediments  more  than  24  cm  deep  in  the  glass  tank  to 
avoid  interference. 

3.  Effects  of  Near  Field  Interference 

In  the  near  field  of  a  transmitter,  the  axial  pressure  exhibits  strong  interference 
effects  (Kinsler.  Frey.  Coppens.  Sanders,  19S2).  Thus,  the  sound  field  is  irregular  and 


does  not  fall  off  smoothly  with  distance  as  it  does  in  the  far  field.  To  avoid 
measurements  in  the  near  field,  an  experiment  was  run  to  find  the  transition  distance 
which  separates  the  near  field  and  far  field.  '1  his  distance  is  20  cm  from  the  transducer 
center.  No  problems  should  occur  as  long  as  the  target  is  more  than  20  cm  from  the 
transducer.  In  the  glass  tank,  we  placed  the  transducer  30  cm  above  sediment.  The 
water 'air  interface  measurement  also  gave  a  check  that  30  cm  distance  is  satisfied. 

Since  a  pulse  moving  in  water  at  14S0  m  s,  takes  405  fts  to  travel  the  30  cm  to  the 
sediment  interface  and  back,  the  30  cm  distance  is  long  enough  so  that  the  transducer, 
which  is  on  about  1SS  fis,  stops  ringing  before  the  first  echo  from  the  sediment  arrives. 

4.  Effects  of  Water/ Air  Interface  Interference 

Ideally,  the  transducer  should  be  far  enough  away  from  all  unwanted  surfaces 
to  avoid  interference.  As  the  transducer  faces  down  and  projects  a  major  part  of  its 
energy  to  the  sediment  surface,  some  energy  is  transmitted  to  the  water  surface  and 
reflected  back  toward  the  sediment.  If  the  transducer  is  close  to  the  water  air  surface, 
the  reflected  energy  will  interfere  with  the  transmitted  energy.  Restricted  by  the  size  of 
the  glass  tank,  we  had  to  locate  the  transducer  close  to  the  water  surface.  To  test  how 
important  this  interference  could  be,  the  water  surface  was  agitated  by  fingers.  The 
waveform  displayed  on  both  analog  and  digital  oscilloscope  was  unchanged  during  the 
ruffling.  It  was  concluded  that  this  source  of  interference  is  negligible  for  this 
experiment. 


41 


5.  Effects  of  Side  Lobe 

The  half  beamwidth  (to  the  first  null)  for  the  F-41  transducer  was  measured  to 
be  10.0°.  For  a  30  cm  distance  from  sediment  to  the  transducer,  the  sediment  surface 
ensonified  by  the  major  lobe  is  a  circle  area  with  radius  5.3  cm.  If  the  distance  is  SO 
cm.  the  radius  will  be  14.1  cm.  For  both  the  glass  tank  and  wood  tank,  if  the 
transducer  is  located  above  the  center  of  the  tank,  the  major  lobe  will  not  produce 
interference  caused  by  the  sides  of  the  tank. 

The  glass  tank  was  too  small  to  avoid  the  interference,  if  any,  caused  by  the 
side  lobe.  To  test  how  much  the  interference  was,  sediments  outside  the  center  30  cm  x 
30  cm  square  area  were  disturbed  and  left  in  a  very  irregular  surface.  Waveforms 
returned  from  sediment  were  measured  before  and  after  the  disturbance.  Data  set  51 
(Figure  4.15)  shows  the  waveform  measured  before  the  sediment  was  disturbed  and 
data  set  52  (Figure  4.16)  shows  the  waveform  after  disturbance.  By  comparing  the  two 
waveforms,  it  is  found  that  the  F-41  trandsucer  gave  negligible  side  lobe  interference  in 
this  experiment. 

6.  Effects  of  Sediment  Inhomogeneity 

Before  the  test  findings  on  side  lobe  effects,  we  rvere  bothered  by  the  drastic 
changes  on  the  waveform  whenever  the  transducer  was  moved  parallel  to  the  smoothed 
aggregate  surface.  It  was  assumed  that  the  waveform  fluctuation  was  caused  by  the 
inhomogeneitics  of  the  aggregate.  To  prove  this,  we  placed  the  transducer  at  a  fixed 
position  in  the  glass  tank  and  repeatedly  ensonified  the  aggregate.  After  each  data  set 
was  acquired,  the  aggregate  was  disturbed  and  smoothed  again.  Ten  sets  of  waveform 
data  were  obtained  and  each  of  them  differed  from  the  others. 

It  was  thought  the  fine  sand  would  be  more  homogeneous  than  aggregate  and 
it  might  produce  more  consistant  waveforms  as  the  transducer  was  moved  horizontally. 
Measurements  made  with  the  fine  sand,  however,  showed  that  the  waveform  still 
fluctuated.  After  each  data  set  was  obtained  from  the  smoothed  fine  sand,  the 
transducer  was  moved  horizontally  to  different  position.  All  the  transducer  positions 
were  at  least  15  cm  from  the  wall  to  avoid  interference  caused  by  reflection  from  the 
sides  of  tank.  Also,  the  positions  were  at  least  a  length  of  tranducer  diameter  from 
each  other  to  avoid  a  waveform  replication.  A  metal  plate  was  set  on  top  of  the  sand 
and  the  transducer  was  moved  horizontally.  As  expected,  the  waveform  on  the 
oscilloscope  was  unchanged  but  totally  different  from  sand. 
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Figure  4.15  Decav  envelope  for  reflection  from  aggrcgatc/watcr 
interlace  in  the  glass  tank  (data  det  5 1)7 
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Figure  4.16  Decav  envelope  for  reflection  from  aggregate, 'water 
interlace  in  the  glass  tank  (data  set  52). 
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CORRELATION  COEFF-  .0400 


E.  PROGRAMMING 


An  HP-S6  desktop  computer  was  adopted  to  collect  and  analyze  the  waveform 
data.  The  computer  program  was  written  to  be  interactive.  This  enables  1)  the 
collection  of  data  sets  with  different  delay  time  intervals  and  2)  the  collecting  and 
processing  of  numerous  data  sets.  The  interactive  programs  allowed  essential 
information  to  be  entered  via  the  key  board,  eliminating  the  need  to  modify  the 
program  whenever  data  sets  were  changed.  The  interactive  prompts  were  useful  for 
people  unfamiliar  with  Basic  programming. 

1.  Program  "THESIS3" 

This  program  was  designed  as  the  initial  step  in  the  data  collection  processing. 
Before  the  selected  waveform  data  displayed  on  the  3091  Nicolet  Digital  Oscilloscope 
were  output  to  the  computer,  it  was  switched  into  "freeze"  or  "store"  mode  from  "live" 
mode.  Through  the  interface  RS-232.  the  computer  received  the  stored  waveform  data 
and  then  processed  them.  After  processing  the  waveform  data  and  before  receiving  the 
next  sampling  data,  the  Nicolet  3091  had  to  be  switched  into  "live"  mode  again.  An 
HP-3421  A  Data  acquisition  and  Control  Unit  was  interfaced  between  the  Nicolet  3091 
and  HP-S6.  The  HP-3421  A  reduced  the  amount  of  manual  operations  in  the  data 
collection  process.  A  five  volts  source  were  input  to  the  HP-3421  A.  Commands  form 
the  HP-S6  and  5  volts  signal  were  used  to  control,  by  closing  and  opening  the  channel 
in  the  HP-3421A,  the  "mode  switching"  of  the  oscilloscope.  To  avoid  errors,  the 
oscilloscope  should  be  warmed  up  enough  while  entering  the  inital  information.  Once 
this  initial  information  was  entered  and  the  oscilloscope  was  warmed  up  enough:  the 
whole  system  runs  automatically. 

Waveform  data  output  from  the  oscilloscope  were  a  stream  of  characters.  To 
obtain  the  voltage  and  time  value,  all  these  characters  needed  to  be  normalized  and 
converted  by  the  computer.  A  program,  listed  in  the  Nicolet  3091  operation  manual, 
was  used  for  data  value  normalization.  A  data  set,  the  waveform  of  the  decay  portion, 
was  obtained  by  squaring  and  averaging  50  samples  of  waveform  data  for  each 
transducer  position.  The  running  average  formula  used  in  the  program  was 

<  V2>n  =  ((n-1) <  V2  >n_!  +Vn2)  n 

where  <  V2  >  is  n  samples  averaged  squared-voltage,  <  V2  >a_j  is  n-1  samples 
averaged  squared-voltage  and  Vn2  is  squared  voltage  for  the  nth  sample.  Note  that 
the  average  window  for  the  waveform  is  100  fis  from  the  beginning  of  the  decay.  The 
above  formula  was  processed  for  each  corresponding  time  bin. 
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During  each  data  set  processing,  some  error  messages,  such  as  STRING  OVF 
(which  stands  for  the  string  variable  overflow),  were  displayed  on  the  HP-S6  screen. 
This  error  messaged  could  be  fatal  to  the  data  set  if  the  data  has  not  been  stored  on 
the  disk.  To  avoid  this  problem,  the  program  was  designed  so  that  after  every  10 
samples  processed  the  updated  data  were  stored  on  the  disk  over  the  previous  updated 
data.  The  program  eliminated  the  problem  of  the  disk  being  full  before  the  final  data 
set  was  stored  on  the  disk.  This  prevention  could  increase  the  processing  time,  but 
compared  to  the  time  wasted  from  the  above  problems,  it  is  worthwhile.  The  final 
result  for  this  program  was  3  data  files  :  one  stored  the  squared  voltage,  one  for  time 
data  and  one  for  the  number  count  of  the  samples. 

2.  Program  "MAXIMUM" 

In  order  to  obtained  the  envelope  of  the  waveform,  this  program  was  designed 
to  pick  up  the  local  maximum  of  the  voltage  data.  The  voltages,  stored  in  the  data  file 
from  "THESIS3",  were  read  and  compared  with  the  neighboring  voltage  data  to  decide 
the  local  maximum.  The  determined  local  maximums  and  corresponding  time  bins 
were  stored  in  seperated  files  and  printed  out.  The  data  files  were  differentiated  by 
number  according  to  the  order  they  were  acquired. 

3.  Program  "PLOT3" 

After  the  local  maximums  were  determined,  the  decay  portion  of  the  waveform 
can  be  plotted  by  program  "PLOT3''.  The  experiments,  for  this  project,  took  place 
mainly  in  the  glass  tank.  A  decay  envelope  measured  in  the  glass  tank  has  an  interval 
between  520  and  619  fis  starting  from  transmitting  time.  Some  exceptions  existed  and 
gave  different  time  frames. 

To  determine  the  decay  constant,  the  voltage  data  were  converted 
logarithmically  and  plotted  vs.  the  time 

V  =  VQ  e-«t 

LN  (  V  )  =  -at  +  LN  (  VQ  ) 

If  the  cursor  on  the  oscilloscope  was  not  centered,  the  resulting  envelope  showed  a 
"zig-zag"  instead  of  a  straight  line  for  a  perfect  exponential  decay. 

The  least-square  method  was  adopted  to  calculate  the  best  fit  of  straight  line 

y  =  ax  +  b 
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where  a,  which  is  similar  to  the  decay  constant,  is  the  slope  and  b  is  the  intercept  on 
the  y-axis.  According  to  Beers  (1953), 


a  =  - 


kI  - 1  xn  I  >n 

kI  w  -  <1  xn  r 


b  =• 


I  *„2 1  yn  -  y  Xny  (xnyn) 

kI  xn2  -  (I  xn  )2 


where  k  is  the  total  number  of  pairs  of  (x,y)  values,  n  is  the  number  between  1  and  k, 
Y  xR  is  the  sum  ot  the  total  x  value,  V  yn  is  the  sum  of  the  total  y  value.  V  xn~  is  the 
sum  of  the  squared  x  value.  X  yn“  is  the  sum  of  the  squared  y  value  and  X  (xnyn)  is 
the  sum  of  the  multiples  of  xn  and  v  . 

Note  that  Y  xn“  and  T  (xn)“  are  not  the  same.  V  x  ~  implies  that  each 
value  of  x  is  squared  and  then  a  sum  is  made  of  these  squares,  while  T  (xn)2  implies 
that  all  of  the  values  of  x  added  together  and  then  this  sum  is  squared.  There  is  an 
analogous  distinction  between  X  (xn'n)  an<^  X  xnS  -  n'  our  case<  tke  ^-coordinate 
is  the  time  value  and  the  y-coordinate  is  the  natural-log  of  the  voltage.  To  check  how 
close  the  observed  and  least-square  estimated  v-value,  compare,  a  correlation 
coefficient  was  calculated  (Chatfield,  19S4) 


(  xn  '  x  )(  -  n  -  y  ) 

r  - - 

(  I  (  xn  •  x  <  -  n  •  >’  y  ) 12 

where  r  is  the  correlation  coefficient,  x  is  the  mean  value  of  x  and  v  is  the  mean  value 
of  y.  Slope  and  correlation  coefficient  were  both  printed  with  the  graph. 

4.  Program  "AVERAGE" 

The  received  echo  waveform  varies  from  position  to  position  because  of  the 
inhomogeneous  sediment.  To  find  the  decay  constant,  we  have  to  average  out  the 
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inhomogeneities  to  obtain  a  mean  value  of  the  waveform.  This  program  performs  a 
running  average  on  the  original  data  sets,  producing  the  final  averaged  waveform  of  the 
sediment. 

The  above  four  programs  are  given  in  Appendix  A. 
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V.  RESULTS 


The  data  sets  for  reflection  from  aggregate  water  interface  are  data  42  to  51. 
Figure  C.l  to  Figure  C.9,  and  Figure  4.15.  Table  3  shows  the  slope  obtained  for  each 
data  set  varied  between  data  sets.  To  get  a  mean  value  for  the  slope,  the  running 
average  method  was  used  to  process  the  data  sets.  The  running  average  method  also 
shows  how  fast  the  average  processing  became  stable.  Data  sets  DAT2  to  DAT10. 
Figure  C.10  to  Figure  C.1S,  are  the  results  for  each  aserage  process.  The  slopes  are 
listed  on  Table  4.  It  is  evident  that  the  slope  becomes  stable  after  S  sets  of  data  are 
averaged.  The  final  a  verged  slope  is  2.13  x  10'“  Np  ps  and  is  quite  different  from  the 
value  reflected  from  the  water  air  interface.  The  difference  between  the  two  slopes 
resulted  from  the  effects  of  volume  reverberation  in  the  aggregate,  and  the  scattering  is 
significant. 


TABLE  3 

ORIGINAL  AGGREGATE  REFLECTION 


DATA 

SLOPE 

CORRELATION 

SET 

( xlO”2 ) 

COEFFICIENT 

DATA42 

-0.  S4 

-0.  5278 

DATA43 

-1.  35 

-0.  5874 

DATA44 

-2.  68 

-0.  7848 

DATA45 

-2.  64 

-0. 8509 

DATA46 

-0.  12 

-0. 1162 

DATA47 

-1.  81 

-0. 3971 

DATA48 

-0.  95 

-0. 3990 

DATA49 

-3.  77 

-0. 8654 

DATA50 

-2.  47 

-0.  9756 

DATA51 

-0.  16 

-0. 0584 

The  data  sets  for  reflection  from  fine-sand  water  interface  are  data  55  to  64, 
Figure  D.l  to  Figure  D.10.  The  slope  of  each  data  set  is  listed  on  Table  5.  Evidently. 
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TABLE  4 

AVERAGED  AGGREGATE  REFLECTION 


DATA 

SLOPE 

CORRELATE 

SET 

( Xl0~2) 

coeeeici: 

DAT2 

-1.  19 

-0.  6051 

DAT3 

-1.  97 

-0.  9720 

DAT  4 

-2.  32 

-0. 9175 

DAT5 

-2.  18 

-0.  9343 

DATS 

-2.  10 

-0.  9914 

DAT7 

-1.  98 

-0. 9952 

DATS 

-2.  16 

-0. 9932 

DAT9 

-2.  20 

-0. 9943 

DAT  10 

-2.  13 

-0. 9958 

the  waveform  changes  very  much  for  different  data  sets.  The  same  running  average 
method  was  used  to  remove  the  inhomogenities.  Data  sets  DAT  12  to  DAT20  were 
produced  from  the  average  process.  They  are  listed  on  Table  6  and  plotted  from 
Figure  D.l  1  to  Figure  D.19.  The  slope  became  stable  after  6  data  sets  were  averaged. 
This  might  suggest  that  the  fine  sand  has  a  greater  homogeneity.  The  averaged  slope 
for  the  fine  sand  is  2.37  x  10"2  Xp  ps  and  is  greater  than  the  aggregate.  However,  the 
slope  difference  between  them  is  not  as  distinguished  when  compared  to  the  slope 
difference  between  the  sediment  and  the  water  air  surface. 

According  to  Equation  3.1,  the  decay  portion  volume  scattering,  represented  in 
voltage,  for  the  aggregate  is 

Va=  K(e-21300t  -  e-37700t) 

and  for  the  line  sand  is 

V  —  K  {  e'237^1  -  e'377OOt  ^ 

The  above  relations  only  existed  without  surface  scattering  or  the  surface  scattering 
was  negligible. 
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TABLE  5 

ORIGINAL  FINE  SAND  REFLECTION 


DATA 

SLOPE 

CORRELATE 

SET 

( xiO”2 ) 

coeffici: 

DATA55 

-1.  51 

-O. 8178 

DATA56 

-3.  03 

-0.  9581 

DATA57 

-2,  50 

-0.  8112 

DATA58 

-0.  75 

-0. 7290 

DATA59 

-1.  00 

-0.  6776 

DATA 60 

-3.  41 

-0.  9957 

DATA 61 

-3.  33 

-0. 9050 

DATA62 

-1.  42 

-0. 4922 

DATA63 

-2.  41 

-0. 9389 

DAT AS 4 

-1.  03 

-0. 5784 

A  simplistic  model  for  sediment  volume  reverberation  was  developed  and  is 
described  in  Appendix  A.  Combining  the  experimental  results  and  the  newly  developed 
model,  the  volume  scattering  coefficient  obtained  was  0.0413  ±  0.007  m°  for  fine  sand 
and  0.0624  ±  0.007  m  for  aggregate.  These  value  are  reasonable  compared  to  the 
difference  between  the  decay  constant  obtained  from  the  experiment.  They  are 
produced  as  a  best  fit  using  the  least  square  method.  The  correlation  coefficients, 
approximately  0.88,  indicate  they  are  quite  a  good  fit.  As  a  first  model,  the  result  is 
promising.  However,  the  discrepancy  between  the  experiment  and  model  might  suggest 
the  model  is  still  very  crude.  Also  mentioned  in  the  Appendix  A  is  that  this  is  only  a 
first  step  to  try  to  obtain  a  theory  to  describe  the  volume  scattering  efi'ect  of  the 
sediment:  it  need  to  be  improved  by  further  effort. 

The  ka  value  for  fine  sand  is  0.12,  which  is  quite  small  compared  to  1.  From 
Equation  2.2,  cloud-of-small-spheres  model,  the  volume  scattering  coefficient  for  the 

'y 

fine  sand  is  0.0082  m°  and  is  not  too  far  from  0.0413  m  .  As  mentioned  in  Chapter 
2,  Equation  2.2  was  derived  under  various  assumptions  and  does  not  apply  to  real 
sediment.  From  Equation  2.3,  single-small-sphere  model,  the  volume  scattering 
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TABLE  6 

AVERAGED  FINE  SAND  REFLECTION 


DATA 

SLOPE 

CORRELATE 

SET 

( xlO“2) 

COEFFICL 

DAT  12 

-1.  84 

-0. 9265 

DAT  13 

-2.  09 

-0. 9051 

DAT  14 

-1.  46 

-0.  9309 

DAT  15 

-1.  47 

-0. 9223 

DAT  16 

-2.  36 

-0.  9912 

DAT  17 

-2.  39 

-0.  9904 

DAT  18 

-2.  39 

-0. 9902 

DAT  19 

-2.  38 

-0. 9903 

DAT20 

-2.  37 

-0.  9916 

coefficient  is  0.34  m'J  which  is  an  order  of  magnitude  different  from  the  the  measured 
value  0.0413  m°.  These  values  suggest  that  the  cloud-of-small-spheres  model  is  closer 
to  the  newly  developed  model  than  the  single-small-sphere  model. 

The  scattering  coefficient  obtained  from  Equation  2.2  is  9456  and  from 
Equation  2.3  is  1134  m°.  Obviously  the  two  values  deviate  considerably  front  the  true 
value.  The  ka  value  is  equal  to  2.1,  which  is  outside  the  Rayleigh  scattering  limit. 

Thus  Equations  2.2  and  2.3  can  not  be  applied  to  the  aggregate. 
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Vi.  CONCLUSION  AND  RECOMMENDATION 


The  following  conclusions  are  possible: 

1.  The  inhomogenities  within  the  sediment  caused  the  reflected  waveform  to 
change  considerably  for  various  transducer  positions.  The  changes  can  be 
removed  by  averaging  over  the  inhomogeneities  and  the  averaging  yields  a 
mean  decay  constant. 

2.  The  giass  tank  works  very  well,  and  therefore  a  bigger  tank  to  contain  the 
sediment  is  not  necessary. 

3.  The  laboratory  procedure  to  measure  the  volume  reverberation,  using 
off-the-shelf  equipment  and  a  high  pass  filter  in  a  water-over-sediment 
system,  is  workable. 

4.  The  basic  program,  written  for  a  HP-86  desktop  computer,  is  sufficient  for 
this  experiment. 

5.  Although  Equations  2.2  and  2.3  can  not  be  applied  to  aggregate  since  it 
has  a  ka  greater  than  1.  the  aggregate  does  have  a  higher  volume 
reverberation  than  fine  sand.  The  aggregate  has  a  larger  size  than  the 
fine  sand;  this  might  suggest  that  the  larger  size  sediment  causes  more 
volume  reverberation. 

6.  From  the  scattering  coefficient  0.0413  ±  0.007  m‘~  for  the  fine  sand  and 
0.0624  ±  0.007  m’“  for  the  aggregate,  the  scattered  energy  received  by 
transducer  is  quite  small  compared  to  the  incident  energy. 

7.  For  fine  sand  sediment,  the  cloud-of-small-spheres  model  works  better 
than  the  single-small-sphere  model. 

8.  The  simplistic  model,  derived  from  the  true  acoustical  pulse  and  described 
in  Appendix  A,  gives  a  crude  formulated  volume  scattering  relationship. 
Combining  with  the  experiment  result,  the  sediment's  volume  scattering 
coefficient  can  be  estimated. 

Further  study  of  the  volume  reverberation  from  ocean  sediment  is  recommended. 
Future  experimentation  should  include  further  variation  of  the  sediment  size  and  of  the 
acoustic  frequency  to  find  a  relationship  combining  the  ka  value  and  volume 
reverberation.  The  developed  model  should  be  improved  to  obtain  a  better  result 
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coincident  with  the  experimental  result.  The  absorption  coelTicient  of  different 
sediment  with  the  frequency  used  should  be  measured  to  compare  the  value  obtained 
from  Urick's  formula  described  in  appendix  A.  Equations  2.2  and  2.3  should  be 
improved  to  obtain  the  scattering  coefficient,  which  can  be  used  to  check  the  above 
model.  Continued  use  of  the  experimental  procedure,  equipment  and  computer 
programs,  written  for  this  project  on  the  HP-S6,  is  highly  recommended. 
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APPENDIX  A 

DERIVATION  OF  SCATTERING  FROM  THE  SEDIMENT 

The  generated  electrical  pulse  can  be  described  by  the  step  unit  function.  II,  as 
follows: 


I(t)  =  IQ  (  H(t)  -  IJ(t-T)  )  (A.l) 

where  I(t)  is  the  intensity  at  time  t,  IQ  is  the  intensity  amplitude,  T  is  the  pulse  length, 

0  t  <  0 

H(t)  = 

1  t  ^  0 

and 

0  t-r  <  0 

II(t-T)  = 

1  t-r  >  0 

Since  the  transducer  is  a  resonant  system  with  damping,  the  square  electrical  pulse  was 
transformed  into  an  acoustic  pulse  with  an  exponential  rise,  an  approximately  flattened 
top  and  exponential  decay.  The  shape  of  the  acoustic  waveform  can  be  described, 
(Sanders  and  Coppens,  19S6),  as 

f(t)  =  (1  -  e'at)  H(t)  +  (1  -  e"ar  )  Il(t-r)  (A. 2) 

where  t  is  the  delay  time  and  a  is  the  decay  constant. 

Let  a  transducer,  located  in  the  water,  send  out  a  sound  pulse  to  a  smoothed 
sediment  surface.  To  assure  the  surface  scattering  is  negligible  and  the  volume 
scattering  is  only  single  scattering,  reflected  pulse  received,  by  the  same  transducer, 
should  have  an  intensity  (Sanders  and  Coppens,  1986) 

I  ~  (J  Tj2  sv  e-Ps  io  fft-2S.'c)d£/(r+S)2  +  Ri  KO/r2  )  /  r2  (A.3) 
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Where  T:  is  the  intensitv  transmission  coefficient.  s,r  is  the  volume  scattering 
coefficient,  P  is  the  absorption  coefficient  in  the  sediment,  £  is  the  affected  sediment 
depth.  Rj  is  the  intensity  reflection  coefficient  and  r  is  the  distance  between  the 
transducer  and  sediment.  From  above  equation,  the  intensity  contributed  by  scattering 
is 

Is  =  Jff  sv  e'^  I0  f(t-2u  c)d£  (r  +  £)2  r2  ( A.4 ) 


and  the  intensity  contributed  by  the  reflection  is 
Ir  =  Rj  10  f(t)  r4 
If  r  >  >  q.  then 

ls  =  (Tj2  sv  I0  r4  )Je-^  Rt-2'  c;dq 


(A.  5) 


(A. 6) 


where 


flt-2^c)  =  (1  -  cMt-2%'c) )  H(t-2^  c)  -  (1  -  )  H(t-t-2£'c)  (A. 7) 


During  the  integration,  the  lower  boundary  of  £  is  zero;  the  upper  boundary  is  ct  2  for 
H(t-2^  c)  and  c(t-r)  2  for  H(t-r-2q  c).  Note  that  only  the  decay  portion,  t  >  t.  was 
interested.  After  the  integration, 

Is  =  A  (  -(e'cPt/2  -  l),p  -  (e*cPt/2*  e‘at)/(-P  +  2a  c)  (A. 8) 

+  (e-cP^‘T)  -  1),  p  +  (c-cP^-t)-'2  .  c'a^T)  );(-p  +  2a/c)) 

where 

A  =  Ti2svI0/r4  (A. 9) 

From  the  experimental  formula  p  =  0.25f  (Urick,  1979)  and  fis  182  kHz,  the  absorption 
coefficient 
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P  =  45.5  dB  m  =  5.2  Np  m, 


cP  2  is  4030  Np  s  for  aggregate  and  is  41S6  Np  s  for  the  fine  sand.  The  intensity  decay 
constant  a  is  twice  the  pressure  decay  constant  5.  From  the  experiment.  6  is  37700 
Np  s  and  it  gives  an  a  equal  to  75400  Np  s;  2a  c  is  97.2  Np  m  for  the  aggregate  and  is 
93.6  Np  m  for  the  fine  sand.  After  plugging  in  the  above  values  for  the  aggregate. 


Equation  A.S  becomes 

Is  =  A(e*4030t  11.6)  -  A(S.3  e‘75400t) 

(A.  10) 

For  the  fine  sand,  Equation  A. 7  becomes 


Is  =  A(e'4,S6t  11.0)  -  A(S.6  e'75400t) 

(AH) 

Jp  Jp 
f-  II 

\  j 

II 

The  pressure  reflection  coefficient  R  is  0.356  for  the  aggregate  (Diaz,  19S6)  and  0.36 
for  the  fine  sand  (Bradshaw,  19S1).  For  the  aggregate, 

Rj  =  0.1270 
T;  =  0.8730 

For  the  fine  sand, 

Rj  =  0.1296 
Tj  =  0.S704 

The  pulse  length  is  SS  ps  and  from  equation  A. 5, 


Ir  =  IQ  Rj  (e*at  -1)  e'75400t  /r4 

(A.  12) 

Using  the  above  values,  the  reflected  intensity  for  the  aggregate  becomes 


Ir  =  96.6  e'75400t  1Q  /r4 

(A.  13) 

and  for  the  fine  sand 


Ir  =  98.5  e'75400t  IQ  / r4 

(A.  14) 
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The  intensity  contributed  by  the  scattering  for  the  aggregate  is 


Is  =  (sv  IQ  r4  )(  e'4030t,  15.2  -  6.3  e'754001  ) 


(A.  15) 


and  for  the  fine  sand. 


ls  =  (sv  I0  r4  )(  e'41S6t  14.5  -  6.5  e'75400t  ) 


(A.  16) 


Combining  the  above  equations,  the  reflected  intensity  from  the  aggregate  received  by 
the  transducer  is: 


I  =  d0  r4)  B 


(A.  17) 


where 


B  =  (  sy  e'44~:(-'*L  15.2  -  6.3  sv  e"'4^  4-  96.6  e'/,?4('' 


400 1 


(A.  IS; 


For  the  fine  sand, 

I  =  ( 1 0 ' r4)  D  (A. 19) 

where 

D  =  (  sv  e‘41S6t,  14.5  -  6.5  sy  e'75400t  +  9S.5  e'75400t  )  (A.20) 

By  inputting  a  series  of  scattering  coefficients  and  computing  the  Log  (B)  and 
Log  (D)  ,  the  theoretical  envelope  was  determined.  A  least-square-fit  slope  was 

obtained  for  each  envelope.  The  scattering  coefficient  for  each  sediment  can  be 
estimated  by  the  value  of  sv  for  which  the  slope  of  the  theoretical  envelope  is  equal  to 
the  slope  of  experiment  data.  For  the  fine  sand,  measured  the  decay  constant  23700 
Np/s  gave  a  scattering  coefficient  of  0.0413  nT3.  For  the  aggregate,  the  scattering 
coefficient  0.0624  in  was  obtained  from  the  measured  decay  constant  21300  Np  s. 

Figure  A.l  and  Figure  A. 2  show  the  envelopes  determined  theoretically  from  the 
above  values  of  sy  for  the  aggregate  and  fine  sand  respectively.  The  theoretical 
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envelope  is  plotted  by  normalizing  the  left  end  of  the  best  fit  slope  line  to  the  same 
value  obtained  from  the  experiment  and  then  applying  this  normalization  to  the  whole 
theoretical  envelope. 

Results  from  the  fine  sand  were  selected  to  find  the  uncertainty  of  scattering 
coefficient  obtained  from  the  above  method.  By  comparing  a  series  of  envelope  plots 
for  different  scattering  coefficient,  the  uncertainty  =  0.00“  was  determined.  Figure  A. 3 
and  A.-  show  the  envelopes  for  scattering  coefficients  =  0.00”  different  from  0.0- 1 3 
m°.  The  slope  changes  as  the  scattering  coefficient  changes.  In  the  beginning  of 
decay,  the  dominant  term  in  equation  A. 20  is  the  refection  term  which  is  independent 
on  s...  and  the  scattering  term  which  denends  on  s,.  is  negligible.  The  theoretical  value 
will  not  change  much  in  the  beginning  of  decay  as  the  scattering  coefficient  changes. 
The  normalization,  described  previously,  still  gave  a  good  approximation  to  the  plot  of 
the  theoretical  envelope  vs.  the  experimental  envelope  for  different  scattering 
coefficients. 

Note  that  the  model  derived  and  described  in  this  appendix  is  crude.  A  later 
impro  .emen.  cn  mo lic i  is 
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Figure  A.l  For  the  aggregate,  the  decay  envelope  obtained  from  the  model 
vs.  the  decay  envelope  obtained  from  the  experiment. 
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TIME  (MICRO-SED 


Figure  A. 2  For  the  fine  sand,  the  decay  envelope  obtained  from  the  model 
vs.  the  decay  envelope  obtained  from  the  experiment. 
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Figure  A. 3  For  the  (Inc  sand,  the  decay  envelope  obtained  from  the  model 
vs.  the  decay  envelope  obtained  from  the  experiment. 
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Figure  A. 4  For  Hie  line  sand,  the  decay  envelope  obtained  from  the  model 
vs.  the  decay  envelope  obtained  from  the  experiment. 
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APPENDIX  B 

COMPUTER  PROGRAM 


i  OOO 
1  0  1  0 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1  090 
1  1 00 
1110 
1120 
1 1 30 
1  140 
1  150 
1  1  60 
1170 
1  1 90 
1181 
1  182 
2000 
20 1 0 
2020 
2021 
2022 

2023 

2024 

2025 

2026 

2027 

2028 

2029 

2030 

203 1 

2032 

2033 

2034 

2035 

2036 

2037 
2040 
2060 
2061 
2062 

2063 

2064 
2080 
2100 
2120 
2140 

2219 

2220 


I  I  I  I  I  I  !  I  I  I  I  I  I  I  I  1  I  I  I  I  I  I  I  I  I  I  t  I  I  I  .  1  I  I  I  I  »  I  I  I  I  I  I  I  !  ■  I  I  I  I  |  (  |  I  I  I  |  I  '  |  I  I  l|  » 


PROGRAM  :  THESIS3 


THIS  PROGRAM  IS  DESIGNED  TO  CONVERT  AND  NORMALIZE  THE 
THE  WAVEFORM  DATA  VALUE  AND  POINT  LOCATION,  TRANSFERED 
FROM  THE  NICOLET  3091  DIGITAL  OSC I LLASCOFE ,  TO  THE 
VOLTAGE  AND  TIME  VALUE.  THE  NICOLET  IS  COTROLED  BY 
THE  HP-3041  DATA  ACQUISITION  AND  CONTROL  UNIT. 

THE  PROGRAM  SQUARE  THE  VOLTAGE  VALUE  AND  THEN  TAKE 
THE  RUNNING  AVERAGE  FOR  FIFTY  SET  OF  READING.  TO 
INSURE  THE  DATA  IS  NOT  L03SED  DURING  COMPUTER  ERROR, 
THE  DATA  IS  BEING  STORE  IN  THE  DISK  AFTER  EVERY  10 
SET  CALCULATION. 


i  ;  i  i  i  i  l  i  t  i  l  i  I  i  i  i  i  i  i  i  i  i  i  t  I  i  i  I  I  i  i  i  i  i  t  t  i  i  i  t  i  j  i  i  i  i  i  i  j  i  •  i  j  i  i  i  ;  i  i  t  l  t  i  i 


DIM  0(1501) , T ( 1 50 1 ) ,3(1501) 

DIM  VD AT AT L 20 ] , VT I ME  S C 20 3 , VCOUNT T C 20 3 
INTEGER  N, I 

i 

!  INPUT  THE  FILE  NAME  TO  STORE  DATA 


DISP  "WHAT  IS  THE  FILE  NAME  FOR  VOLTAGE  DATA  ?" 

INPUT  VDATAT 

DISP  "WHAT  IS  THE  FILE  NAME  FOR  TIME  DATA  ?" 

INPUT  VT I  MET 

DISP  "WHAT  IS  THE  FILE  NAME  OF  THE  COUNT  DATA  ?" 

INPUT  VCOUNT T 

DISP  "YOUR  INTEREST  TIME  INTERVAL  IS  FROM  C START  3  TO  CENDT]  " 
DISP  "WHAT  IS  THE  STARTING  TIME  OF  MEASURING  IN  MICRO-SEC)  ?" 
INPUT  START 

DISP  "WHAT  IS  THE  ENDING  TIME  OF  MEASURING  (  IN  MICRO-SEC)  ?" 
INPUT  ENDT 


!  CHECKING  THE  INTERESTED  TIME  INTERVAL  FROM  CRT  SCREEN 

i 

LET  N= 1 

i 

!  IF  THIS  IS  NOT  THE  FIRST  SET  OF  DATA,  THE  COMPUTER  IS  GOING 
!  TO  READ  THE  DISK  FOR  PREVIOUS  DAT  INFORMATION,  AND  CARRY  ON 
!  THE  JOB 

i 

DISP  "IS  THIS  YOUR  FIRST  DATA  SET  OUT  OF  THE  50  DATA  SET?" 

DISP  "PRESS  O  IF  YES  :  PRESS  1  IF  NO" 

INPUT  ANS 

IF  ANS= 1  THEN  GOSUB  READ_DATA 

i 

!  FOLLOWING  IS  THE  COMPUTER  RECEIVE  THE  DATA  FROM  3091  AND 
!  THEN  CALCULATE  THE  VOLTAGE  AND  TIME  VALUE. 


2240 

2260 

2280 

2290 


DIM  I  $-  C  5  3  ,  DT  C  20000 1  ,  NT  C  40 1 
CONTROL  10,3  :  15 

CONTROL  10,4  ;  58 
CONTROL  10,2  :  1 
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2300 

2301 

2302 

2303 

2304 

2305 

2306 

2307 
2300 
2310 
231  1 

2312 

2313 

2314 

2315 

2316 

2320 

2321 


OUTFUT  70?  ; " OF  N4 " 

i 

«  ABOVE  STATEMENT  IS  LET  THE  CRT  TO  '’LIVE1'  CONDITION 

i 

i 

!  FOLLOWING  ONE  IS  TO  LET  THE  CRT  IN  MSTOREM  CONDITION 

i 

t 

OUTFUT  70?  ; "CLS5" 

WAIT  1000 

I 

!  FRESS  BOTTOM  ONLY  USED  WITHOUT  THE  HP-3421  UN  I T 
'  DISP  "FRESS  3091  BUTTON" 

i 

'  LET  THE  DATA  TRASMIT  TO  THE  CONFUTE?  THROUGH?  RS-232 

i 

OUTFUT  709  ;"CLS2n 


ABOVE  STATEMENT  LET  THE  WAVEFORM  INFORMATION  BE  TRANSFERED 
FROM  CRT  TO  COMPUTER  THROUGH  RS-232  INTERFACE 


2340  ENTER  10  :  It 
2360  ENTER  10  :  Dr 
2300  ENTER  10  :  Nr 
2390  1 

2400  HO=VAL  (NtC 16,203) 

2420  V 1  =  ( VAL  ( N  t  C  2 1 , 253 ) -5)  *10'  (VAL  (M t C 26 . 30 3 ) - 1 2) 

24  4  0  HI = (VAL  ( NT  C  3  1 , 35  3  )  -5 )  *  10'x  (VAL  (NT  C  36 , 40  3  )  -  1  2) 

2441  * 

2443  OUTFUT  709  : M0FN5M 

2445  OUTFUT  70?  s "0FN2M 

2446  •  LET  THE  CRT  IN  LIVE  MODE 

2447  OUTFUT  70?  ;"CLS4" 

2460  J  =  1 

2400  K=5 

2500  F  R I NT  "N=",N 
2520  ■  FOR  1=0  TO  ENDT 

2525  IF  P- START  OR  I ; ENDT  THEN  2780 

2530  ! 

2531  !  NORMALIZE  THE  VOLTAGE  VALUE  AND  TIME  FOR  THE  INTEREST  INTERVAL 

2532  ! 

2540  D(I)=VAL  ( D t C J , K 3 ) *V 1 

2560  T ( I ) = ( I —HO ) 1 

2500  D< I )=D( I >  2 

2600  IF  M= 1  THEN  GOTO  2640 

2620  S ( I ) = ( (N- 1)*S(I)+D(I))/N 

2621  1 

2640  IF  N= 1  THEN  S(I)=D(I) 

2660  IF  I =ST ART  THEM  2760  ELSE  2600 
2600  IF  I =ST ART  +  20  THEN  2760  ELSE  2700 

2700  IF  I=START  *-30  THEN  2760  ELSE  2720 

2720  IF  I =ST ART  +  50  THEN  2760  ELSE  2740 

2740  IF  I =ENDT  THEN  2760  ELSE  2700 

2760  DISP  T ( I ) ; " SEC  D ( I ) j "VOLT  "  ;  S  ( I  ) 

2770  FRINT  T(I):"SEC  " : D ( I ) ; " VOLT  " ; S  ( I ) 

2771  { 

2780  J-J+5 

2000  K=K  +  5 
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2920  ME XT  I 
2840  N=MJM 
2950  ■ 

2051  1  LET  THE  DAT A  BEING  STORE  AFTER  EVERY  10  SET  OF  CALCULATION 

2852  ' 

2860  IF  N  MOD  10=1  THEM  2965  ELSE  2240 
2965  IF  N= 1 1  THEN  2970  ELSE  2980 

2970  CREATE  VDATAT,200,B 

2971  CREATE  VCOUNTf, 1 ,8 

2972  CREATE  VTIME 1,200,8 
2980  GOSUB  WRITE_DISK 

2990  IF  N=5 1  THEN  3000  ELSE  2240 
3000  END 

3305  f 

3306  ! 

3307  ! 

3308  !  FOLLOW  I NG  SUBROUTINE  IS  USED  TO  WRITE  THE  DATA  TO  THE  DISK 

3309  ! 

3310  ! 

3500  WRITEJ3ISK: 

3510  PURGE  7DAT A  t' 

3520  CREATE  VDATA* , 200 , 8 
3540  ASSIGNS  1  TO  VDATAf 

3550  PURGE  VCOUMT r 
3560  CREATE  VCOUMT*, 1,8 
3580  ASSIGN#  2  TO  VCOUMT t 

3590  PURGE  VT I  MET 

3600  CREATE  VT I MET , 200 , 8 
3620  ASSIGN#  3  TO  VT I ME * 

3640  PRINT#  2  :  M 

3660  FOR  I  =START  TO  ENDT 

3680  F  R I  NT#  3  :  T(I) 

3700  PRINT#  1  :  S ( I  ) 

3720  NEXT  I 

3740  ASSIGN#  1  TO  * 

3760  ASSIGN#  2  TO  * 

3780  ASSIGN#  3  TO  » 

3900  RETURN 

3900  ! 

3910  ! 

3920  «  FOLLOWING  SUBROUTINE  IS  USFD  TO  READ  THE  DATA  FROM  FREVIOUS 

3930  '  WRITTEN  D I  SI 

3940  ! 

3950  ' 

4000  READ_DAT A: 

4020  ASSIGN#  2  TO  VC0UNT1 

4040  ASSIGN#  1  TO  VDATAT 

4060  READ#  2  ;  N 

4061  PRINT  N 

4080  FOR  I =S  T AR  T  TO  ENDT 
4  100  READ#  1  s  S  < I > 

4120  IF  I =ST ART  THEN  4220  ELSE  4  140 
4140  IF  I =ST  AR T  +20  THEM  4220  ELSE  4  160 

4  160  IF  I -ST ART  +  30  THEN  4220  ELSE  4  180 

4180  IF  I =START  +  50  THEN  4220  ELSE  4200 

4200  IF  I =ENDT  THEN  4220  ELSE  4240 
4220  PRINT  S(I) 

4240  NEXT  I 

4260  ASSIGN#  1  TO  * 

4280  ASSIGN#  2  TO  * 

4380  RETURN 
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300 

801 

802 

803 

804 

805 

806 

807 

808 

809 

810 

811 

812 

813 

814 

815 

316 

880 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

991 

992 

993 

994 

996 

1 OOO 

1 020 

1040 

1060 

1080 

1  1 00 

1  120 

1140 

1  141 

1144 

1146 

1148 

1  150 

1152 

1  154 

1  1 60 

1  165 

1  170 

1  180 

1 200 

1220 

1  240 

1260 

1  280 

1290 

1291 


t  i  ;  i  i  *  I  i  i  i  i  t  i  i  i  i  i  i  i  t  i  i  i  i  i  i  t  i  t  i  i  i  i  i  i  i  i  i  i  <  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  t  i  i  i  i  * 


PROGRAM  :  MAXIMUM 


THIS  PROGRAM  IS  DESIGNED  TO  READ  THE  SQUARED 
WAVEFORM  VOLTAGE  DATA  AND  TIME  DATA  STORED 
IN  THE  DISK  AND  THEN  CHOOSE  THE  LOCAL  MAXIMUM 
(WAVEFORM  ENSEMBLE).  THE  CHOOSED  MAXIMUMS  WILL 
BE  TAKEN  SUARE-ROOT  AND  STORE  THE  FINAL  RESULT 
ON  THE  DISK 


i  i  i  t  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  »  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i 


DIM  VDATA.SC  20]  ,VTIME*C20]  .  VNUMBER  *  [  20  ]  ,  VT  I  ME  _MAX*  C  20  ]  ,  VMAX  I MUMS*  C  20 1 
DISP  "WHAT  IS  THE  NAME  FGF:  VOLTAGE  DATA  FILE  ?" 

INFUT  VDATA* 

DISP  "WHAT  IS  THE  NAME  FOR  TIME  DATA  FILE 
INPUT  VT I  ME  £ 

DISP  "WHAT  IS  THE  NAME  FOR  FILE  TO  STORE  THE  NUMBERS  NO.  ?" 

INPUT  VNUMBER* 

DISP  "WHAT  IS  THE  NAME  FOR  FILE  TO  STORE  THE  T I ME  _MAXS  DATA  ?" 

INPUT  VT  I  ME_MAX-S 

DISP  "WHAT  IS  THE  NAME  FOR  FILE  TO  STORE  THE  MAXIMUMS*?" 

INPUT  VMAX I  MUMS* 

DISP  "WHAT  IS  THE  STARTING  TIME  OF  YOUR  DATA  COLLECING  ?" 

INPUT  START 

DISP  "WHAT  IS  THE  ENDING  TIME  OF  THE  DATA  COLLECTING  ?" 

INFUT  ENDT 


ASSIGN#  1  TO  VDATA* 
ASSIGN#  3  TO  VT I ME  S 
CREATE  VNUMBER*, 1,8 
ASSIGN#  4  TO  VNUMBER* 
CREATE  VTIMEJ1AX*, 100,8 
ASSIGN#  5  TO  VT  I  ME_MAX-S 
CREATE  VMAX I MUMS*, 100,8 
ASSIGN#  6  TO  VMAX I MUMS* 


INTEGER  PRINTTER 

DISP  "WHAT  IS  THE  PRINTER  DEVICE  ADDRESS  MO.  ? 
'  INPUT  THE  PRINTER  NO.  SCREEN  IS  1 
INPUT  PRINTTER 
PRINTER  IS  PRINTTER 


DIM  SC1501) , T ( 1 50 1 ) , TIM (100) ,MAXS(100) 

PRINT  USING  1170  ;  "TIME  (SEC "VOLT " 

IMAGE  25X,  UA,  12X,4A 
INTEGER  I,J 
FOR  I -START  TO  ENDT 
READ#  3  ;  T  (  I > 

READ#  1  ;  S(I> 

NEXT  I 
LET  J  =0 

t 

'  IF  THE  FIRST  DATA  BEING  READ  IS  GREAT  THAN  THE  SECOND  ONE, 
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1292 

1 193 

1 300 

1320 

1340 

1360 

1380 

1  400 

1410 

1420 

1440 

1460 

1480 

1500 

1520 

1540 

1  560 

1 5B0 

1 600 

1610 

1620 

1640 

1660 

16  BO 

1 700 

1720 

1730 

1740 

1760 

1  BOO 

1820 

1840 

1  BBO 

1  900 


!  THE  FIRST  ONE  IS  TREATED  AS  A  LOCAL  MAXIMUM. 

i 

IF  S(START>  '=  S  (ST ART  + 1 )  THEN  1320  ELSE  1500 
T I M ( J ) =  T ( START  > 

MAXS (J) =SQR  (S (START) ) 

PRINT#  5  ;  TIM(J) 

PRINT#  6  ;  MAXS(J) 

PRINT  USING  1410  ;  " J=" ; J , T IM < J ) , MAXS ( J ) 

IMAGE  12X , 3 A , DDD ,  7X , D . DDDDDD , 1 IX ,D. DDDDDDDD 
!  PRINT  T I M ( J ) 

!  PRINT  MAXS(J) 

J=J  +  1 


FOR  I =ST ART  + 1  TO  ENDT-1 

IF  S(I).  =  S(I-l)  AND  S  (  I  >  .:•=  S(I  +  1)  THEN  1540  ELSE  1720 
T  I M  ( J  )  =  T  (  I  > 

MAXS (J) =SOR  (S ( I ) ) 

i 

PRINT  USING  1610  ;  " J-" , J , TIM ( J ) ,MAXS(J) 

IMAGE  12X ,3A,DDD.7X ,D. DDDDDD, 1 1 X , D . DDDDDDDD 
!  PRINT  TIM ( J ) 

1  PRINT  MAXS ( J  > 


PRINT# 

cr 

wJ  « 

TIM(J) 

PRINT# 

6  : 

MAXS ( J) 

J=J  +  i 

NEXT  I 

J=J-i 

PRINT# 

4  : 

J 

ASSIGN# 

4 

TO 

* 

ASSIGN# 

5 

TO 

* 

ASSIGN# 

6 

TO 

* 

ASSIGN# 

1 

TO 

* 

ASSIGN# 

3 

TO 

END 
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900  1  • | i | ■! i | i i i i i i  j  i i i i i i i i  >  i i i i i i i i i i i i  >  i i  »  i i i i i i i i i i i i i i i  j  j i  <  ; i 

901  1  !  .  ! 

903  !  !  ' 

Q04  !  !  PROGRAM  :  PLOT 3  ! 

905  .  (  I 

906  !  '  ) 

907  !  !  THIS  PROGRAM  IS  DESIGNED  TO  PLOT  THE  ENVELOP  ! 

900  !  !  OF  WAVEFORM  BY  READING  DATA  FROM  DATA  DISK.  ! 

909  '  !  AFTER  READING,  VOLTAGE  DATA  IS  CONVERTED  TO  LOG  ! 

Q10  !  !  VALUE,  THEN  BE  PLOTTED.  AFTER  PLOTING,  A  LEAST-  ! 

911  !  !  SQUARE  METHOD  IS  USED  TO  FIND  THE  BEST  FIT  STRAIGHT  ! 

912  !  ■  LINE  FOR  THE  LOG-ENSEMBLE.  PLOT  THE  LINE,  PRINT  ' 

913  !  !  OUT  THE  SLOPE  OF  THE  LINE  AND  THE  CORRELATION  ' 

914  1  !  COEFFICIENT.  ,  ! 

915  !  !  ! 

916  '  1  ! 

917  |  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  t  1  1  1  1  1  1  1  1  t  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  *  1  1  1 

910  ! 

919  ' 

1000  CLEAR 
1020  GCLEAP 

1022  INTEGER  PLOTER 

1024  DISP  "WHAT  IS  YOUR  POLTTER  DEVICE  ADDRESS  NO.  ?“ 

1025  !  INPUT  THE  PLOTTE  NO.  SCREEN  IS  1 

1026  INPUT  PLOTER 
1028  PLOTTER  IS  PLOTER 

1040  FRAME 

1041  ! 

1042  !  THIS  PLOT  SIZE  IS  TO  FIT  THE  REQUIREMENT  FOR  THE  NFS  THESIS 

1043  ! 

1060  LOCATE  30,125,30,85 

1000  CSIZE  3 

1100  MOVE  50,20 

1120  LABEL  “TIME  (MICRO-SEC)" 

1140  PEN  UP 
1160  MOVE  10,40 
1100  DEG 
1200  LDIR  90 

1220  LABEL  "LN <<< VOLTAGE  2>)  1 /2) " 

1240  PEN  UP 

1241  ! 

1242  '  DEFAULT  SET-UP  IS  TO  FIX  THE  PLOTING  COORDINATION  AS 

1243  !  X-AXIS  IS  FROM  500  TO  670  (MICRO-SEC),  Y-AXIS  IS  FROM 

1244  !  -5  TO  2 

1245  ! 

1250  DISP  “DO  YOU  WANT  TO  USE  DEFAULT  SET-UP  ?  YES-0 “ 

1251  INPUT  BULL 

1252  IF  BULL=0  THEN  1253  ELSE  1260 

1253  SCALE  500,670,-5,2 

1255  LAXES  10, .5,500,-5,2,2 

1256  XMIN=500 

1257  GOTO  1710 

1259  ! 

1260  DISP  “ENTER  THE  XMIN  OF  SCALE" 

1200  INPUT  XMIN 

1300  DISP  “ENTER  XMAX  OF  SCALE'1 
1320  INPUT  XMAX 

1340  DISP  "ENTER  THE  YMIN  OF  SCALE" 

1360  INPUT  YMIN 
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1380 

1  400 

1420 

1440 

1460 

1480 

1  500 

1 520 

1540 

1560 

1 580 

1600 

1620 

1640 

1660 

1661 

1680 

1700 

1710 

1720 

1721 

1740 

1741 

1742 

1744 

1746 

1748 

1749 

1750 

1751 

1752 

1759 

1760 

1780 

1781 

1782 

1790 

1800 

1802 

1803 

1 804 

1810 

1820 

1B21 

1822 

1823 

1825 

1826 

1828 

1831 

1832 

1833 

1834 

1  835 

1837 

1839 

1840 

1920 

1930 

1940 


DISP  "ENTER  THE  YMAX  OF  SCALE 
INPUT  YMAX 

SCALE  XM  IN  , XMAX , YM IN , YMAX 
FXD  0,2 


DISP  "ENTER 
INPUT  XT 
DISP  "ENTER 
INPUT  YT 
DISP  "ENTER 
INPUT  XI 
DISP  "ENTER 
INPUT  Y I 
DISP  "ENTER 
INPUT  XMC 
DISP  "ENTER 


THE  X-TICKING  SPACE" 
THE  Y-TICKING  SPACE" 
THE  X  INTERSECTION" 
THE  /  INTER  SECTION" 
THE  X-MAJOR  COUNT" 
THE  Y-MAJOR  COUNT" 


INPUT  YMC 

LAXE3  XT , YT , X  I  , Y I , XMC , YMC 
GRAPH 
PEN  1 


DIM  TIM ( 100) ,MAXS ( 100) 

DIM  TIMEMAX*C203 , MAX  I MUM* £  20 3 

DISP  "WHAT  IS  YOUR  FILE  NAME  FOP  TIMEMAX*  ?" 

INPUT  TIMEMAX* 

DISP  "WHAT  IS  YOUR  FILE  NAME  FOR  MAXIMUM*  ?" 

INPUT  MAXIMUM* 

PRINT  "THE  DATA  FILES  APE " , T I MEMAX* , " AND " , MAX  I MUM* 

DISP  "WHAT  IS  THE  HIGHEST  ORDER  OF  ARRAY  NUMBER?" 

!  THIS  GIVES  A  CHOOICE  TO  CHOOSE  THE  END  FOINT  OF  PLOT 
INPUT  JK 

i 

ASSIGN#  5  TO  TIMEMAX* 

ASSIGN#  6  TO  MAXIMUM* 

DISP  "PRESS  O  IF  THE  FIRST  LOCAL  MAXIMUM  IS  BEING  SAVED  ?" 

INPUT  OK 

FOR  I«0  TO  JK 

READ#  5  :  TIM(I) 

1  CONVERT  THE  TIME  UNIT  FROM  SEC  TO  MICRO^SEC 

TIM(I)=TIM(I)-*1 000000 
READ#  6  ;  MAXS(I) 

MAXS ( I ) =LOG  (MAXS(I) ) 

!  CONVERT  THE  VOLTAGE  TO  NATURE  LOG  VALUE 
NEXT  I 
K-JK+ 1 

IF  0K=0  THEM  GOTO  1835 
FOR  1=1  TO  JK 
TIM ( I  —  1 ) =T I M ( I ) 

MAXS ( I  —  1 ) =MAXS ( I ) 

NEXT  I 
K=JK 

1  DISP  "WHICH  CHARACTER  YOU  PREFER  TO  USING  IN  PLOTTING  THE  DATA? 
!  INPUT  C* 

MOVE  TIM (0) , MAXS (0) 

FOR  1=0  TD  K—l 
DRAW  TIM ( I )  , M A X S ( I ) 

!  LABEL  C* 

NEXT  I 
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1960 

1961 

1962 
I960 
2000 
2020 
2900 
2910 
2920 
2922 
2924 
2926 
2928 
2930 
3000 
3040 
3060 
3080 
3 1 00 
3120 
3140 
3160 
3  1  50 
3200 
3220 
3240 
3260 
3280 
3300 
3320 
3340 

3350 

335 1 


3360 

3370 

337 1 

3372 
3380 
3400 

3420 

3421 

3422 

3423 

3440 

3441 

3442 

3443 

3460 

3461 

3480 

3481 
3500 

3520 

3521 
3540 

3560 

3561 

3562 
3580 


FEN  UP 

GOSUB  LEASTSQR 
ASSIGN**  5  TO  4 
ASSIGN#  e>  TO  •* 
END 


‘  FOLLOWING  SUBROUTINE  IS  DESIGNED  TO  APPLY  THE  LEAST 

!  SQUARE  METHOD  IN  THE  ABOVE  DATA  AND  THE  PLOT  THE 

!  BEST  FIT  STRAIGHT  LINE,  PRINT  OUT  THE  SLOPE  AND 

!  THE  CORRELATION  COEFFICIENT  ON  THE  PLOT. 

!  FORMULA  USED  P LEASE 

t 

LEASTSCR: 

REAL  Y  ,  X  Y  ,  X  1 , X 2 , TAVE , YAVE , VARX , VARY , CQVXY 
LET  Y=0 

!  Y  IS  THE  SUMMATION  OF  Y<I>  VALUE 

LET  X Y=0 

!  XY  IS  THE  SUMMATION  OF  X < I ) *Y ( I > 

LET  X 1=0 

!  XI  IS  THE  SUMMATION  OF  X(I) 

LET  X2-0 

!  X2  IS  THE  SUMMATION  OF  X ( I )  2 

FOR  1=0  TO  K-l 
Y=Y+MAXS ( I ) 

XY=XY+TIM ( I ) *MAXS ( I ) 

X 1  =  X 1 +T I M  (  I ) 

X2=X2+T I M  (1)2 
NEXT  I 

B=  (X24Y-X  1*XY)  /  (K*X2-X  1  "'2) 


•  A  IS  THE  SLOPE  OF  THE  FITTED  STRAIGHT  LIGHT 

A= ( K  4  X  Y - X 1 4  Y ) / (K*X2-X1  2) 


!  B  IS  THE  INTERCEPT  TO  Y  AXIS 

PRINT  "THE  SLOPE  A=" ;A 
PRINT  "THE  Y- INTERCEPT  B=";B 
YAVE* (MAXS  (0) +  MAXS (K-l )  ) /2 

i 

!  YAVE  IS  THE  AVERAGE  VALUE  FOR  THE  LOG  VOLTAGE 

i 

T  AVE= (TIM ( O ) +T I M ( K— 1 )  ) /2 

i 

!  TAVE  IS  THE  AVERAGE  VALUE  FOR  THE  TIME  VARIABLE 

i 

COVX Y=0 

i 

!  COVX Y  IS  THE  COVARIANCE  OF  TIME  AND  VOLTAGE  (LN) 

i 

VARX=0 

!  VARX  IS  THE  VARIANCE  OF  TIME  VARIABLE 

i 

VARY=0 

!  VARY  IS  THE  VARIANCE  OF  VOLTAGE (LN)  VARIABLE 

DIM  YL(IOO) 

FOR  1=0  TO  K-l 
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7600 

3620 

3640 

3660 

3600 

3700 

3720 

3740 

3760 

3761 

3730 

3800 

3820 

3840 

3860 

3880 

3900 

7910 

3920 

3940 

3960 

3980 

3990 

3991 

3994 

3996 

3998 

3999 

4000 

4002 

4004 

4006 

4020 

4021 

4040 

4060 

4061 

4062 

4080 


C0VXY=C0VXY+  (TI!i<  I  )  -TAVE )  ■»  (MAX  3  (  I)  -YAVE) 

i 

VARX=VARX+ (TIM ( I ) -TAVE)  '2 

i 

VARY=VARY+  (MAXS  <  I )  -YAVE)  '-2 


YL  ( I >=A*TIM  < I ) +8 

!  YL  IS  THE  VALUE  FITTED  ON  THE  STRAIGHT  LINE 


NEXT  I 
PEN  2 

MOVE  TIM (0) , YL <0> 

FOR  1=0  TO  K-l 
DRAW  T I M ( I ) ,YL ( I ) 

NEXT  I 

CORCOE=COVX Y/SQR  (VARX+VARY) 

i 

!  CORCOE  IS  THE  CORRELATION  COEFFICIENT 

i 

FRINT  "CORRELATION  COEFFICIENT  IS  CORCOE 

PEN  l 

DEG 

LDIR  0 

i 

!  FOLLOWING  PRINT  OUT  THE  SLOPE  AND  COORELATION  COEFFICIENT 
!  ON  THE  PLOT 

i 

MOVE  XMIN+80,1 

!  LABLE  ALWAYS  STARTS  FROM  THE  XMIN+80  IN  X  COORDINATES 

i 

LABEL  USING  4021  ;  "SLOPE  «" ,A 
IMAGE  7A , DDD . DDDD 
MOVE  XMIN+80,2 

LABEL  USING  4061  ;  "CORRELATION  COEFF= ", CORCOE 
IMAGE  19A , DDD. DDDD 
PEN  UP 
RETURN 
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BOO 

810 
S20 
8  TO 
S40 
850 
860 
870 
880 
890 


■  i  i  i  i  i  i  i  i  i  i  i  <  i  i  i  i  >•<  i  i  i  i  i  i  t  i  i  i  t  i  i  i  i  i  i  i  i  i  i  i  i  t  i  i  i  i 


PROGRAM  :  AVERAGE  1 

THIS  PROGRAM  15  DESIGNED  TO  AVERAGE  ! 

THE  DATA  SETS  TO  GET  A  AVERAGED  ! 

WAVEFORM  DATA.  ! 

i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  <  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i 


1000  DIM  VDATTC253 , VDATATC253 , VDAT1TC253 
1020  DIM  VI ( 1500) ,V2< 1500) ,V3( 1500) 

1040  D IBP  "WHAT  15  THE  ORIGINAL  DATA  FILE  ?  DAT#" 

1060  INPUT  VDATT 

1080  DI5P  "WHAT  15  THE  NEW  DATA  FILE  ?  DATA#" 

1100  INPUT  VDATAT 

1120  DISP  "WHAT  15  THE  AVERAGED  DATA  FILE  ?  DAT 1 # " 

1140  INPUT  VDAT1T 

1160  ASSIGN#  1  TO  VDATT 

1180  ASSIGN#  2  TO  VDAT  AT 

1200  CREATE  VDAT IT , 8 , 1 10 
1220  ASSIGN#  3  TO  VD4T1T 

1 230  ! 

1231  !  N  IS  THE  NUMBER  HOW  MANY  DATA  SET  BEING  AVERGED 

1232  ! 

1240  DISP  "WHAT  IS  YOUR  AVERAGE  No.  N  ?" 

1260  INPUT  N 

1261  ! 

1262  !  AVERAGE  ONLY  APPLY  TO  THE  SAME  ARRAY  NO. 

1263  ! 

1270  FOR  1=522  TO  621 

1271  READ#  1  :  VI (I) 

1272  NEXT  I 

1274  FOR  1=522  TO  621 

1275  READ#  2  ?  V2 ( I ) 

1276  NEXT  I 

1280  FOR  1=522  TO  621 

1300  ! 

1301  !  RUNNING  AVERAGE  IS  ADOPTED 

1302  ' 

1340  V3(I)=( (N-1>*V1 (I>+V2(I> ) /N 

1360  PRINT#  3  ;  V3 ( I > 

1371  ! 

1372  !  PRINT  OUT  FOLLOWING  VALUE  AS  A  CHECKING  PURPOSE 

1373  ! 

1380  IF  1=522  THEN  1500  ELSE  1400 

1400  IF  1=553  THEN  1500  ELSE  1420 

1420  IF  1=600  THEN  1500  ELSE  1440 

1440  IF  1=621  THEN  1500  ELSE  1520 

1500  PRINT  VI ( I > , V2 ( I ) , V3 ( I ) 

1520  NEXT  I 

1540  END 
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APPENDIX  C 

DECAY  ENVELOPE  GRAPHS  FOR  AGGREGATE 

This  appendix  includes  the  graphs  of  the  decay  envelopes  for  reflection  from 
aggregate  water  interface  within  the  glass  tank.  Figure  C.l  to  C.9  show  the  original 
decay  envelopes  obtained  from  different  transducer  positions.  The  decay  envelopes  are 
also  described  as  following  data  sets: 

data  set  42 
data  set  43 
data  set  44 
data  set  45 
data  set  46 
data  set  47 
data  set  48 
data  set  49 
data  set  50 

Figures  C.10  to  C.18  arc  the  decay  envelopes  resulted  from  the  running  average 
method.  They  are  also  described  as  following  data  sets: 

data  set  DAT2 
data  set  DAT 3 
data  set  DAT4 
data  set  DAT5 
data  set  DAT6 
data  set  DAT7 
data  set  DATS 
data  set  DAT9 
data  set  DAT  10 
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Figure  C.  1  Decav  envelope  for  reflection  from  aggregate, 'water  interface, 
transducer  was  30  cm  above  surface  (data  set  42). 


77 


LN(«V0LTAGE‘2>)  *1/2) 


2 

m 


r> 

X) 

o 

i 

cn 

rn 

n 


i  i  i  i  *  i 


I'igure  C.2  Decay  envelope  for  reflection  from  aeeregate, 'water  interface, 
transducer  was  30  cm  above  surface  (data  set  43). 
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TIME  (MICRO-SEC) 
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CORRELATION  COEFF=  -.7848 


TIME  (MICRO-SEC) 


Figure  C.4  Dccav  envelope  for  rellection  from  agercgatc/watcr  interface, 
transducer  was  30  cm  above  surface  (data  set  45). 
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CORRELATION  C0EFF=*  -.8509 


TIME  (MICRO-SEC) 
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Figure  C.5  Dccav  envelope  for  reflection  from  aggregate, 'water  interface, 
transducer  was  30  cm  above  surface  (data  set  46). 
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Figure  C.6  Decay  envelope  for  reflection  from  aggregate/water  interface, 
transducer  was  30  cm  above  surface  (data  set  47). 
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Figure  C.7  Decav  envelope  for  reflection  from  aggregate/water  interface, 
transducer  was  30  cm  above  surface  (data  set  48). 
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Figure  C.S  Decay  envelope  for  reflection  from  ageregate/water  interface, 
transducer  was  30  cm  above  surface  (data  set  49). 
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Figure  C.9  Dccav  envelope  for  rcllcciion  from  aggregate/water  interface, 
transducer  was  30  cm  above  surface  (data  set  50). 
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f  igure  C.  10  Decay  envelope  resulted  from  the  average  of  data  set  42  and  43 
['or  the  aggregate/water  surface  (data  set  DA *12). 
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Figure  C.l  1  Decov  envelope  resulted  from  the  average  from  data  set  42  to  44 
lor  the  aggregate,  water  surface  (data  set  DA  13). 
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Figure  C.  1 2  Decay  envelope  resulted  from  the  average  from  data  set42  to  45 
for  the  aggregate/water  surface  (data  set  DA'1 4). 
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Figure  C.13  Decay  envelope  resulted  from  the  average  from  data  set  42  to  46 
for  the  aggregate/water  surface  (data  set  DAT5). 
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Figure  C.14 


Decay  envelope  resulted  from  the  average  from  data  set  42  to  47 
lor  the  aggregate/water  surface  (data  set  DA'1 6). 
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Figure  C.  15  Dccav  envelope  resulted  from  the  average  from  data  set  42  to  48 
for  tire  aggrcgatc/watcr  surface  (data  set  DA  17). 
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Figure  C.16  Dccav  envelope  resulted  from  the  average  from  data  set  42  to  49 
(or  the  aggregate/water  surface  (data  set  DATS). 
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Figure  C.17  Decav  envelope  resulted  from  the  average  from  data  set  42  to  50 
for  the  aggregate/water  surface  (data  set  DA  19). 
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Figure  C.18  Decay  envelope  resulted  from  the  average  from  data  set  42  to  51 
for  tile  aggregate/water  surface  (data  set  DA'1'10). 
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APPENDIX  D 

DECAY  ENVELOPE  GRAPHS  FOR  FINE  SAND 


This  appendix  includes  the  graphs  of  the  decay  envelopes  for  reflection  from  fine- 
sand,  water  interface  within  the  glass  tank.  Figure  D.l  to  D.10  shew  the  original  decay 
envelopes  obtained  from  different  transducer  positions.  The  decay  envelopes  are  also 
described  as  following  data  sets: 

data  set  55 
data  set  56 
data  set  57 
data  set  58 
data  set  59 
data  set  60 
data  set  61 
data  set  62 
data  set  63 
data  set  64 

Figures  D.ll  to  D.19  are  the  decay  envelopes  resulted  from  the  running  average 
method.  They  are  also  described  as  following  data  sets: 

data  set  DAT  12 
data  set  DAT  13 
data  set  DAT  14 
data  set  DAT  15 
data  set  DAT  16 
data  set  DAT  17 
data  set  DAT  18 
data  set  DAT  19 
data  set  DAT20 
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Figure  D.l  Decay  envelope  for  reflection  from  fine-sand/water  interface, 
transducer  was  30  cm  above  surface  (data  set  55). 
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Figure  D.2  Decav  envelope  for  rellcction  from  fine-sand, 'water  interface, 
transducer  was  30  cm  above  surface  (data  set  56). 


97 


CORRELATION  C0EFF-  -.9581 


TIME  (MICRO-SEC) 


LN  ( (<V0LTAGE'2>)  *  1/2) 


I  I 

LD  ^ 

5C0  p— -i  T" 

I 

520  j- 

i 

\ 

I- 

i 

540  j- 

t 

i 

560 J- 
i 

[ 

| 

580  j- 

1- 

i 

600  [- 

T 

620  [- 

r 

640  1- 


I 

(jl) 


I 

ro 

•~r’ 


T 


ro 

n 


i 

ao 

%  -* 
i—* 

ro 


660  r 


Figure  D.3  Dccav  envelope  for  reflection  from  fine-sand/water  interface, 
transducer  was  30  cm  above  surface  (data  set  57). 
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Figurc  D.4  Decay  envelope  for  rcllcction  from  fine-sand, 'water  interface, 
transducer  was  30  cm  above  surface  (data  set  58). 
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Figure  D.5  Decay  envelope  for  reflection  from  fine-sand/watcr  interface, 
transducer  was  30  cm  above  surface  (data  set  59). 
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Figure  D.6  Decay  envelope  for  reflection  from  fine-sand/water  interface, 
transducer  was  30  cm  above  surface  (data  set  60). 
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Figure  D.7  Decay  envelope  lor  reflection  from  line-sand, 'water  interlace, 
transducer  was  30  cm  above  surface  (data  set  (51). 
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Figure  D.8  Dccav  envelope  for  reflection  from  fine-sand  water  interface, 
transducer  was  30  cm  above  surface  (data  set  62). 
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Figure  D.9  Decav  envelope  for  reflection  from  fine-sand, 'water  interface, 
transducer  was  30  cm  above  surface  (data  set  63). 
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Figure  D.10  Decay  envelope  for  rcllection  from  fine-sand/water  interface, 
transducer  was  30  cm  above  surface  (data  set  64). 
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Figure  D.l  1  Decay  envelope  resulted  from  the  average  of  data  set  55  and  56 
lor  the  fine-sand/water  surface  (data  set  DAT  12). 
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Figure  D.12  Decay  envelope  resulted  from  the  average  from  data  set  55  to  57 
for  tile  fine-sand/water  surface  (data  set  DA  1  13). 


107 


TIME  CMJCRO-SEO 


Figure  D.13  Decav  envelope  resulted  from  the  average  from  data  set  55  to  58 
for  the  (ine-sand, 'water  surface  (data  set  DAT  14). 
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Figure  D.14  Decav  envelope  resulted  from  the  average  from  da 
for  the  line-sand/ water  surface  (data  set  DAT  15). 
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Figure  D.  16  Dccav  envelope  resulted  from  the  average  from  data  set  55  to  61 
lor  the  line-sand/watcr  surface  (data  set  DAT  17). 
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Figure  D.17  Decay  envelope  resulted  from  the  average  from  data  set  55  to  62 
lor  tile  line-sand, 'water  surface  (data  set  DA  I  18). 
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Figure  D.18  Dccav  envelope  resulted  from  the  average  from  data  set  35  to  63 
For  tlic  fme-sand/water  surface  (data  set  DAT  19). 
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I  igure  D.19  Decay  envelope  icsultcd  from  1  he  average  from  data  set  55  to  6d 
lor  t lie  aggrcgatc/watcr  suiiacc  (data  set  DA  I  2<»). 
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